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ABSTRACT 



The models of the non-linear optics in which solitons were appeared are 
considered. These models are of paramount importance in studies of non- 
linear wave phenomena. The classical examples of phenomena of this kind 
are the self-focusing, self-induced transparency, and parametric interaction 
of three waves. At the present time there are a number of the theories based 
on completely integrable systems of equations, which are both generations 
of the original known models and new ones. The modified Korteweg-de 
Vries equation, the non- linear Schrodinger equation, the derivative non- 
linear Schrodinger equation, Sine-Gordon equation, the reduced Maxwell- 
Bloch equation, Hirota equation, the principal chiral field equations, and the 
equations of massive Thirring model are gradually putting together a list of 
soliton equations, which are usually to be found in non-linear optics theory. 

Keywords: solitons, self-induced transparency, non-linear fibres, wave in- 
teraction PACS number: 42.65. 



1 Introduction 



A non-linear wave is one of the fundamental objects of nature. They are in- 
herent to aerodynamics and hydrodynamics, solid state physics and plasma 
physics, optics and field theory, chemistry reaction kinetics and population 
dynamics, nuclear physics and gravity. All non-linear waves can be divided 
into two parts: dispersive waves and dissipative ones. The history of in- 
vestigation of these waves has been lasting about two centuries. In 1834 
J.S. Russell discovered the extraordinary type of waves without the disper- 
sive broadening. In 1965 N.J. Zabusky and M. D. Kruskal found that the 
Korteweg-de Vries equation has solutions of the solitary wave. This solitary 
wave demonstrates the particle-like properties, i.e., stability under propaga- 
tion and the elastic interaction under collision to one another. For these 
reasons these solitary waves were named solitons. In succeeding years there 
has been a great deal of progress in understanding of soliton nature. Now 
solitons have become the primary components in many important problems 
of non-linear wave dynamics. It should be noted that non-linear optics is the 
field, where all soliton features are exhibited to a great extent. 

The self-focusing of light beams and self-induced transparency phenomenon 
set the bright example of an important role of solitons in non-linear optics. 

The non-linear wave propagation in non-resonant medium also occupies 
an important place in non-linear optics. As c consequence of modulation 
instability, continuous wave transforms to number of single pulses. These 
pulses can under certain conditions evaluate into solitons. Nowadays the 
optical solitons propagation in fibres has attracted considerable attention. 

The other classical problem of non-linear optics is the parametric inter- 
action of waves. The harmonic generation, stimulated by both Raman and 
Brillouin scattering, parametric amplifications, sum-frequency mixing, and 
four-wave mixing have been the subject of many investigations in this field. 
The three-wave interaction process is of particular interest because it gives 
us a new instance of the application of the soliton occurrence in non-linear 
system without dispersion. 

There are cases of the Hamiltonian systems when the canonical transfor- 
mation makes the equations of motion trivially integrable after the conver- 
sion to new variables. Then they say that Hamiltonian system admits the 
action-angle variables. If the action-angle variables exist, then the Hamilto- 
nian system is completely integrable. There is powerful tool to investigate the 
integrable systems. It is Inverse Scattering Transform (1ST) introduced by 
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Gardner, Greene, Kruskal and Miura (1967). The implementation of 1ST in 
non-linear optics results in new examples of the optical soliton phenomena. 
Moreover, this method leads to new analytical approaches for near integrable 
systems. 



2 Classical examples 

Let us consider three phenomena mentioned above for that theories have 
been based on completely integrable equations. 



2.1 Self-induced transparency 

A self-induced transparency (SIT) phenomenon consists in the propagation 
of a powerful ultra-short pulse (USP) of light through a resonance medium 
without the distortion and energy loss of this pulse [JTJ El ^, 0|. This phe- 
nomenon is characterised by the continuous absorption and re-emission of 
electromagnetic radiation by resonant atoms of medium in such a manner 
that steady-state optical pulse propagates. In the ideal case the energy dis- 
sipation of the USP is invisible and the state of the resonant medium is not 
varying. In this means the medium is transparent. The group velocity of 
such steady-state pulse, called 27r-pulse or soliton of SIT, is less than the 
phase speed of light in a medium. The group velocity depends on a 27r-pulse 
duration: the shorter is the duration, the higher is its speed [0, ^ |], [5| . When 
two pulses of the different velocities spread in the medium, the second pulse 
may overtake the first and a collision will take place. After the collision, the 
solitons keep their shape and velocity (but in general all other parameters of 
solitons may alter). This fundamental property of the SIT solitons has been 
studied many times both theoretically and experimentally @, ||, [7j . 

From the mathematical point of view this property is a consequence of the 
complete integrability of the reduced Maxwell-Bloch equations, describing 
the SIT in the two-level media with non-degenerated levels M - |TJ|. The 
27r-pulses answer the single- soliton solutions of these equations, and the 
process of "collision" reflects the evolution of the double-soliton solution - its 
asymptotically transformation into a pair of solitons under certain conditions 
(see, for example, || and Jl4| - fll5l). 



The simplest theory describing the self-induced transparency phenomenon 
was developed by McCall and Hahn [fLl ||. In general, the theory of the in- 
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teraction of radiation with an ensemble of two-level atoms is based on the 
Bloch equations for atoms and the Maxwell equations for the classical elec- 
tromagnetic field. In an isotropic dielectric the set of Maxwell equations 
reduced to one equation for the electric field E = IE. For a plane wave with 
constant polarisation vector I one can obtain the following system of total 
Maxwell-Bloch (MB) equations 



d 2 E 1 d 2 E _ Ann A d I d 2 r l ' 
Ik? ~ &~cW ~ c 2 \~W 



(1) 



dr x dr 2 2d dr 3 2d 

^ = -^ ^T = UJari + T Er ^ ^r = -T Er ^ (2) 

where d is the projection of a matrix element of the dipole operator on the 
direction of I, ha is the concentration of resonant atoms. It should be noted 
that the components of Bloch vector n, r 2 , and r 3 depend on the atomic 
resonance frequency u a . Hereafter the angular brackets represent summation 
over all the atoms characterised by the frequency u a . 

The Bloch equations contain products of the field E and the polarisa- 
tions r 2 and r 3 responsible for interference between the opposite propagated 
waves. It has been shown JT6, [T7||, however, that if the density of resonant 



atoms is small enough to make the parameter Aim Ad 2 /hu a less than unity, 
interference may be neglected. It was found that for a typical value of d ~ 
1 Debye, uo a ~ 10 15 s~ l and ha <C 10 23 cm -3 one may not take into account 
the backward wave generation by a forward running pulse. Thus, the MB 
equations convert into the simple reduced Maxwell-Bloch (RMB) system of 
equations: 

dE IdE 2ixn A d I drA 

!h + c~dt~ c~\~^t/ , () 

dn dr 2 2d dr 3 2d 

-8^ = -^ ^ = UJan + T Er3 ' ~~dt = % ' (4) 
It should be emphasised that in both the MB and RMB equations symbol E 
denotes the real value of the electric field strength. However, every so often 
the electromagnetic wave can be represented as a quasi-monochromatic one 



E(z, t) = 2A(z, t) cos[A;o.2 — uJot + <p(z, t)] = £(z, t) exp[i(koz — u>ot)], (5) 
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where ui is the radiation frequency, k is the wave number, and the real 
envelope A(z,t) and the phase p(z,t) are slowly varying functions of z and 
t. This is an approximation which means that the envelope and the phase 
obey the inequalities 



\dA/dt\ < u \A\, \dA/dz\ < k \A\, \d<p/dt\ < u \(p\, \d(p/dz\ < k \p\. (6) 

Besides, the envelope amplitudes are usually so weak that Rabi frequency 
(max\dA/h\ ) turns out to be much less than the resonance transition fre- 
quency. The resulting system of equations can be represented as 

dq ldq , (dtp ldtp\ , 

d- z + cm = - a{p) > q [-o- z + -cTt = a(g) ' (7) 



dQ (. M n dP (. dp\ „ n OR n fa . 

where Au = {uj a — uq), a' = litoj^nA d 2 /hc, q = dA/h is a normalised slowly 
varying pulse envelope and P, Q and R are connected with initial Bloch 
vector components by the relations 



r 1 = P(z, t) sin[& z — u t + p(z, t)] + Q(z, t) cos[fc z — u Q t + p(z, t)], 
r 3 = -R(z,t) 

To describe the SIT phenomenon McCall and Hahn used the system (0) 
and (pj). Furthermore, if we confine the analysis to situations when the input 
optical pulse does not carry any phase modulation, i.e. dp/dz = dp/dt = 
at z = and the form factor of the inhomogeneous line is a symmetrical 
function of frequency detuning Auj, then equations (|7|) and (H) yield dp/dz = 
dp/dt = at any z and t. In this case equations (|7|) and @ reduce to the 
system of SIT equations 
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If an absorption line is homogeneously broadened and exact resonance con- 
dition holds, then the SIT equations reduce to the well-known Sine-Gordon 
equation 



d 2 u 



(11) 



drdC, 



+ sinw = 0, 



where r — (t — z/c), ( = a'z, and q = du/dr. 

These equations (i.e., (|3| - £|), (0 - §), and (?? - ??)) may be represented 
as the condition of the integrability of some linear equations that provides to 
the solution of these equations be the 1ST method. If we assume, that before 
the arrival of an ultra-short pulse all two-level atoms are in the ground state 
and after the passing of the USP all atoms recover in the initial states, then 
boundary conditions are 



Both RMB-equations and SIT-equations with these boundary conditions can 
be solved by the 1ST method in a regular way || - ||15|1 . In general case one 
obtain the iV-soliton solution of these equations. This solution represents L\ 
single solitons and L2 breathers (so that N = L\ + 2L2 ). Breather (or bion - 
soliton-antisoliton bounded state) is an extremely stable solitary wave with 
internal oscillations. It has the same collision stability as ordinary solitons 
in both the bion-bion and the bion-soliton collisions. It is worth to note that 
the Sine-Gordon equation has the same soliton and breather solution. 

It should be note that one-soliton solution of the RMB equations corre- 
sponds the USP without carrier wave and represents a unipolar spike of elec- 
tromagnetic radiation. Sometimes these pulses are named the video pulses. 
Two-soliton solution of the RMB-equations describes the collision of two 
video pulses in the same fashion as it was done for two soliton solution of 
SIT equations by McCall-Hahn. However, the breather solution of the RMB- 
equations can be used to obtain the generation of McCall-Hahn 27r-pulses. 
In || it was shown that the RMB breather is a real solution in the form of 
localised pulse with internal oscillations. Hence, this is an exact analogue of 
the McCall-Hahn 07r-pulse. If the frequency of internal oscillations increases, 
then the envelope of RMB breather can be described by the soliton solution 
of the SIT equations to a high accuracy. Thus admittedly, the 27r-pulse of 
McCall-Hahn is the limiting case of the 07r-pulse of the RMB equations. 



lim R(r,(;u a ) 



1, lim r-1,2 (t, C;^a) = 0. 
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2.2 Optical solitons in fibres 

It is well known that the rate of the information transfer by means of fibre 
optical communication systems (FOCS) with the mode pulse-code modula- 
tion is limited mainly by the effect of dispersion of group velocities. The 
influence of this effect can be suppressed efficiently and ideally can be com- 
pletely excluded if one uses sufficiently powerful pulses of light. Because of 
the non-linear effect of self-influence such pulses in the certain conditions are 
transformed in solitons and their propagation in FOCS does not accompany 
by dispersive spreading. 

It is important to emphasise that soliton does not exist in real communi- 
cation systems in the true sense of the word. The influence of group velocities 
dispersion of the high orders, optical loss and some other effects break a dy- 
namic balance between the non-linear compression of pulse and its dispersion 
broadening. As a result of these effects the optical pulse suffers an envelope 
distortion and damping. But the distance, passed by the soliton in a fibre 
considerably exceeds one, which weak pulse should be passed in the linear 
regime of propagation. The experiments made with powerful optical pulses 
confirm this. We can consider the soliton as a good approach for the real 
non-linear pulses in a fibre under certain conditions. 

The suggestion to use optical solitons for the information transfer along 



the fibre was made in the works [18], [19| and it was demonstrated [p0| , [21], |22 
Later on the equation describing the optical pulse propagation in a fibre with 
account of only the second-order group-velocities dispersion was received in 
works [T2^, and was investigated in |25|, In the soliton theory it was 
known as the non-linear Schrddinger (NLS) equation. 

Let q is the normalised slowly varying complex envelope of the optical 
pulse defined by the following expression 

E(x, y, z, t) = A q(z, t)^(x, y) exp[i((3 z - u Q t)], 

where j3 is the propagation constant depending on the frequency of the 
carrier wave uq, ty(x,y) is a mode function that determines the transverse 
distribution of the electric field over the fibre cross- sect ion. Slowly varying 
envelope of optical pulse is governed by the follow equation |23|, [24| (NLS 
equation): 

^ + s + ^ = °- (12) 
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Here £ = z/Ljj, r = {t — z/v g )t~Q are normalised independent variables of co- 
ordinate and time, accordingly, t p0 is a pulse duration at z = and u fl is the 
group velocity of an optical pulse. The term in (JT2j) with the second derivative 
with respect to r describes the pulse dispersion broadening (s = —1 for the 
normal dispersion and s — +1 for the anomalous dispersion). Lp is the 
dispersion length L D = 4/3 t 2 (\d 2 (3 / dui 2 ]) -1 . The effective refractive index 
n e ff is defined as (3{ui) = (uj/c)n e ff. The third term in is responsible 
for the self-action effect. Coefficient /x is equal to the ratio of the dispersion 
Ld length to the Kerr length Lk, where Lk = c 2 P (2iiujlAl\xR^f /I) -1 ■ Here 
XReff * s the effective non-linear susceptibility responsible for the Kerr effect. 

The complete integrability of the NLS was established in the classical 
works [^, |27|| . Besides the multi-soliton solutions the NLS equations 
has the new type of the solutions named multiple-pole solutions ||. These 
solutions correspond to the multiple points of a discrete spectrum of the 
Zakharov-Shabat problem in 1ST method. It is difficult to realise these solu- 
tions in practice because any small stir of the initial conditions will remove 
degeneration in the multiple points of a discrete spectrum. It means that 
only the exclusive initial optical pulses can be transformed into multiple- 
pole optical solitons. Different properties of the NLS equations are described 
in the excellent books 



2.3 Interaction of three waves 

One of the broadest classes of phenomena in non-linear optics is the transfor- 
mation of the frequency of an electromagnetic radiation propagating in the 
non- linear medium. Harmonics generation of the fundamental wave (pump), 
sum-frequency and difference-frequency mixing are classified among these 
phenomena pi| . Under sufficiently high intensity of a pump the polarisa- 
tion of a medium is not a linear function of the electric field strength of the 
wave. If the frequencies of an electromagnetic field are not in resonance with 
atomic transition frequencies, one can use a standard perturbation theory to 
reveal this dependency. So we can expand polarisation P in a power series of 
electrical field strength. The coefficients of this series are tensors of the nth 
rank x i named as non-linear susceptibility, describe different processes of 
the electromagnetic waves interaction. The non-linear effects, described by 
the nth rank tensors of non-linear susceptibility, are often interpreted as the 
interaction of the (n + 1) waves. 

Let the non-linear characteristics of a medium be described by non-linear 



2S] 



8 



susceptibility of the second order ■ They call it quadratic non-linear 
media. Let the waves with the carrier frequencies u)\ and u> 2 propagate along 
z axis. As the polarisation is the non-linear (quadratic) function of the 
electrical field strengths, the waves with carrier frequencies uj = uj\ ± UJ2, 
uj = 2uj\ and uj = 2uj2 appear in such medium. These waves, in their turn, 
can cause generation of new waves with the frequencies uj = 2u)\ ± 0J2, uj = 
uji ± 2uj 2 , and so on. But in a dispersive medium all these processes are 
not equally efficient. There is a condition of phase matching, which selects 
a certain type of interaction of three waves, leaving all other unaffected. 
Sometimes such phase matching takes place for the waves propagating in the 
same direction. In this case they say about collinear parametric interaction. 
In this case the distance where the interaction of waves occurs can be made 
sufficiently long and, consequently, the effective frequency transformation 
will take place. On the contrary, when the phase matching is achievable only 
for the waves propagating in different directions, their interaction occurs 
only in the field of overlapping of the wave beams. Non-collinear parametric 
interaction is worth to draw special attention as the number of interesting 
results concerning the integrability of three-wave mixing equations have been 
found. 

Let Si, 8 2 and £3 are the slowly varying envelopes of the interacting pulses. 
Let us consider the situation when only the collinear propagating wave with 
sum-frequency or difference frequency is generated. In the slowly varying 
envelopes and phases approximation, the system of equations describing the 
interaction of the three waves can be written in the following in unified form 



' d_ Id 

dz v 2 dt 



d 1 d\ 

— + — — \q 2 = icrqlql exp(+zA/^) (13) 



dz + V 3 di J q3 = - i(T Qi£ e M+ iA kz) 

where a = v /7T7i7i, and j n = Anuj n x i2) (^i,^2)/cn(uj n ), n = 1,2,3. In 
these equations we used Ak = k 3 — (ki + k 2 ), for a sum-frequency mixing 
process uj = uji+uj 2 , and Ak = k 3 — (ki — k 2 ), for difference-frequency mixing 
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lo = uj\ — uj2- Here v n are the group velocities of the corresponding wave. 
The effects of group-velocity dispersion are neglected. The uniformization 
variables (71,2,3 related with 81,82 and 83 as follows 8\ = v/7i?i, 8 2 = v/72?2 
and 83 = y/^q^, , for the case of sum- frequency mixing, and 8\ = — v/7i?3, 
£1 = \/72<Z2j ^3 = y/^hli f° r the case of difference-frequency mixing (the 
indices at the velocities of the first and the third waves should be changed in 
this case). 

We would like to point out that the resonance Raman scattering under 
condition of weak variation of the energy levels population of a medium and 
the scattering of optical waves by an acoustical wave can be considered from 
one position as the specific realisations of the three-wave interaction. In both 
cases the systems of equations, describing these processes, can be transformed 
into one universal system (|13|). Furthermore, one can demonstrate that under 
certain conditions the reduced Maxwell-Bloch equations appear here, which 
describe the propagation of the ultra-short pulse in a resonance medium. 
Due to this property the Raman scattering can be analysed in terms of 1ST 
method. 

The equations for 3-wave interaction (|13"D permit both the infinite num- 
ber of conservation laws and the Backlund transformation. They can be 
presented as the Hamiltonian equations by employing the r-matrix. It was 
shown that these equations pass the Painleve test and that there is a class of 
self-similar solutions, expressed in terms of Painleve transcendents P-V and 
P-VI. 

It is remarkable that the system of equations, describing the parametric 



interaction of three waves, can be solved by 1ST method for a 3D case |37|, |38[ , 
whereas the most of soliton equations are one-dimensional ones. 

It should be emphasised that the 3-wave interaction gives an example of 
the dispersionless propagation of the non-linear waves. They often say that 
soliton is the result of the compensation of dispersion broadening and non- 
linear compression of the wave packet. The 3-wave interaction just demon- 
strates the narrowness of this statement. Due to the absence of the phase and 
group velocities dispersion the solitons in this process do not detach from the 
non-soliton part of the solution (which is often named radiation). It seems 
very difficult to study the process analytically, so one has to confine the in- 
vestigation to some particular solutions. The non-collinear second harmonic 
generation provides an example of a specific case of the 3-wave interaction, 
where an exact solution was found without employing the 1ST method. The 
solution obtained explicitly illustrates the non-separability of the soliton and 
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non-soliton parts of the solution of these three wave interaction equations. 



3 New examples of the integrable systems 

3.1 Generalisation of the self- induced transparency the- 
ory 

The development of the SIT theory is characterised by its going beyond the 
framework of two-level approximation and by the spectral composition of 
the USP field becoming more complex. The latter means that the reso- 
nant medium interacts with radiation containing several carrier frequencies. 
Besides, such factors as the direct interaction between resonant atoms, non- 
linear properties of the dielectric doped with resonant atoms and polarisation 
of the electromagnetic field should be taken into account. 

Let the optical pulse propagate along axis z and the electric field strength 
be presented in the form E = £ (t, z) exp(— iuot + ik§z) + c.c. The carrier 
frequency is in resonance with the frequency uj 2 \ = (W 2 ~ W\)/h of an atomic 
transition j 2 — > ji between energy levels W 2 and W\ degenerated over the 
projections m and I of the total angular momenta ji and j 2 . In a general case 
the evolution of the envelope of the USP and states of the resonant medium 
are described by the system of equations for which the exact solution is not 
known in the case of arbitrary valuesji and j 2 . But the certain choice of 
transitions ji = <-> j 2 = 1, ji = 1 <-> j 2 = 1, and j x = 1/2 <-> j 2 = 
1/2 makes this system of equations exactly integrable. Its solution can be 
obtained by the inverse scattering transform method, as it was shown in 
HO, Generalised reduced self-induced transparency equations (GSIT 



equations), for transitions ji — <-> j 2 — 1, ji — 1 j 2 — 1 5 can be 
represented in the homogeneous form 



dq 



E Pa (Pj 



=1,2 



|- - tAcot p0 ) P\ a) =^ qMf - Qj N^ , (14) 
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If we consider transition ji = — > j 2 = 1, then in (|T3D we assign g,,- = 
dSHp /h, (3 1 = 1 and /3 2 = 0. For transition = 1 — > j 2 = we have 
<7j = dSHpo/h, j3 l = and /3 2 = 1 . Finally, for transition ji = 1 — > j 2 = 1 we 
have qj = jd£H p0 /h^/2, /3 1 — f3 2 — 1/2- Everywhere here the sub-index and 
upper index take the values j = ±1. The slowly varying envelopes of matrix 
elements of the density matrix p are determined as P„- = (j2,Q\p\ji, j) , 

p} 2) = (h,-j\p\ji,o), = (j 2 ,o\p\j 2 ,o), iv (2) = -(ji,o\P\h,o), M f = 

(2) 

\32, j\p\ji,l) , = (j2, —l\p\ji, —j). As we are well aware, GSIT equations 
are the zero-curvature condition in the vector expansion 1ST method for the 
AKNS hierarchy. The spectral problem with such a sort of 1ST was first 



reported by Manakov |42| to describe the self- focusing of the polarised light 



beams. The expressions for iV-soliton solutions, breathers and Backlund 



transformation have been found in 43 



There are versions of the GSIT equations describing USP propagation 
in three-level medium. In the simplest case of this model the resonance 
levels have V and A configurations. It was determined H3], ESI that, if the 



oscillator forces for every transition in a V or A configuration are equal, then a 
two-frequency pulse (characterised by two different frequencies of the carrier 
wave) is able to propagate in such a medium without the envelope distortion. 
An ultra-short pulse of the kind was called simulton |[46|| . The simultons 
are single-soliton solutions of the above GSIT equations. At the same time 
solutions occur which are responsible for the propagation and collisions of 
simultons. The oscillating simultons (colour breathers) are the two-frequency 
generalisations of 07r-pulse by McCall-Hahn. It is worth noting at this point 
that a simulton is generally unstable with respect to transformation into one- 
frequency 27r-pulse and may remain as two-frequency pulse only for a special 
choice of resonance level populations. 

It should be mentioned that consideration of the polarised USP propaga- 
tion ion three-level medium leads to matrix variants of the GSIT- equations. 



For more details, see Ref. [47, 48]. The matrix expansion 1ST method for 



the AKNS hierarchy (or matrix Manakov spectral problem) there occurs. 

3.2 Femtosecond optical solitons in fibres 

In order to describe the non-linear phenomena associated with femtosecond 
optical solitons in non-linear fibres, the higher-order non-linear Schrddinger 



(HNLS) equation has been proposed f|9l (50J (see also 01). Consider the 
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following generalisation of the equation (|T2] ) 



■ 9q d 2 q . l2 ./ d 3 q . . 2 dq d\q\ 2 \ n 

Z 9C + ^ + ^ 9 + * + 9r + ) = ( } 

The parameter rj 3 corresponds to the third-order group-velocity dispersion, 
parameters fi 2 and /x 3 represent the two inertial contributions to the non- 
linear polarisation, i.e., Raman self- scattering and self-steeping formation. 
If rj 3 = 0, /i 2 =/i 3 =l, and /i=0, then the HNLS flT5|) reduced into derivative 
non-linear Schrddinger (DNLS) equation 

s^ + i\q\ 2 q) = 0. (16) 




This equation is completely integrable one p2 . There are soliton and multi- 
soliton solutions, which can be obtained in the framework of 1ST method. 
Notice that the condition fi=0 is not essential for reduction of the HNLS into 
the integrable equation, i.e., resulting modified DNLS equation is completely 
integrable too. 

If rj 3 = 1, /i 2 = ±6, /i 3 = 0, then the equation (]T5|) is reduced to the Hirota 
equations |BB[ , 



9( 1 , Q2< i I I 1 2 I n , 1<7 N 

% d( + s d^ + g + ^ lgl ^ = ' ( } 

which represents other example of completely integrable equation. 

In general, equation (|15"D may not be completely integrable. However, if 
we suppose that r/ 3 — e, // 2 = 6e, // 3 = 3e:, and // = 2s, then the HNLS 
equation can be reduced to the other one - Sasa-Satsuma equation. In this 
case in terms of new variables £ = r — s 2 (/3e, q(r, £) = w(£, C) exp(isr/3e + 
2is 3 (/27e 2 ) the equation ([15]) takes the form 

d 3 u , l9 du d\u\ 2 \ , , 

^ + %| 2 ^ + 3»^)=0. (18) 

The equation ( |T8| ) has been considered in (||54|), where it was shown that 
this equation could be solvable by means of 1ST method. 
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3.3 Vector optical solitons 

In the general case the vector soliton is the soliton solution of the non-linear 
system of the evolution equations, which one can be presented as the one- 
dimension array. For example, the vector soliton of the NLS equation is the 
solution of the following vector equation 

*f + ^0 + ^tf= ' (19) 

where q = {qi, q%, qu}- Hereafter, this equation will referred to as the 
v-NLS equation. The vector index of this soliton can be resulted from the 
different physical origin. 

We have only one example where the v-NLS equation is the completely 
integrable one. This equation is embedded into AKNS hierarchy that allows 
to exploit the 1ST method to find the soliton solution. The suitable spectral 



problem was found by Manakov [fl2]| . 




As the simplest vector generation of the DNLS equation one can write 
the standard form of the v-DNLS equation 

- ie(q?)?j = 0. (20) 

Other example of the vector non-linear waves arises under consideration 
ultra-short optical pulse propagation in birefringent fibres with higher order 
effects like the third order dispersion of group velocities, Kerr dispersion, and 
stimulated Raman scattering. There is two- component generalisation of the 
Sasa-Satsuma equation (ft 



+ W^ + %T)S + 3^]=0. (21) 



' d^q dq d(qq* 

In the three-component generalisation of this equation was considered 



too. It was shown that these equations have the zero-curvature representation 
and can be solved by 1ST method. However, the exact soliton solutions more 
simply to obtain using the Darboux- Backlund transformation. 
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3.4 Extremely short pulse propagation in non-resonant 
medium 



The recent progress in the field of generation of femtosecond pulses has made 
it necessary to revise theoretical models of their propagation in a non-linear 
dispersive medium. It is interesting to find a method to describe an USP 
evolution without the use of slowly varying envelope approximation. The 
simplest way to do this is to combine the wave equation for electromagnetic 
field with the equations specifying the changes in the state of the medium. 

In the case of resonant medium we have the resonance transition frequency 
as a scale parameter for time. When pulse duration t v obeys the inequality 
tpU) a ^> 1, the slowly varying envelope approximation is adequate to describe 
the pulse propagation. On the contrary, if t p oj a C 1 we can use at least the 
unidirectional propagation approximation. The ratio e = ujR/uj ai where ojr is 
the Rabi frequency (i.e., ujr = dmax.\E\/h), provides a new parameter. Let 
the USP amplitude (i.e., E) be of such magnitude that the Rabi frequency 
ujr is small compared with the minimum of atomic transition frequency. 
That means e = uon/uj a is a small parameter. Then we can attempt to 
solve the Bloch equations (^) approximately and thus obtain an approximate 
equations of the USP electric field strength without the assumption of slowly 
varying envelope p6| . The Maxwell equations under unidirectional wave 



approximation can be reduced into 

BE IdE 2irn A d /dP 1 \ 

dz + cdt ~ c \ dt J ' 1 ] 

where polarisation of the ensemble of two-level atoms at the third order of e 
is 



Substitution of the expression for polarisation into equation fl22|) in terms of 
new variables r = \b\z, ( = t — z/V, u(t, £) = —(a/6bY^ 2 E(z,t) yields the 
modified Korteweg-de Vries equation (mKdV) 

^ +ftI Va? =0 ' (24) 

Here parameters a = ( 2A-nnA\d\ A / cU^uj a \ b = (4imA\d\ 2 /ch 3 uj a ) were in- 
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troduced. The expression V^ 1 = c _1 [l + (4imA\d\ 2 /hu a )} defines the re- 
normalised velocity of the USP propagation. As it is known [|57j] this equation 
is a completely integrable, and its solutions can be found by the 1ST method 
|3T| , |32| , |33|j . There are another non- linear equations describing extremely 
short pulse propagation H5S| , however, they have not soliton solutions. 



3.5 SIT in a Kerr-type non-linear medium 

There are intensive investigations of the non-linear pulse propagation in the 
optical fibres |3(J . However, practically all materials used for fibre fabri- 
cation contain impurities that contribute to the absorption spectrum of the 
fibres. The losses due to the resonant absorption decrease if the frequency 
of the carrier wave is located within the window of transparency of glass fi- 
bre. Another means to decrease losses is to make pulse duration shorter than 
the characteristic relaxation times of the resonant states or in other words 
to make optical pulses ultra-short. In this case the well-known self- induced 
transparency phenomenon can be expected to arise. 

It is known that the non-linear Schrodinger equation, which is used to 
describe optical solitons in non-linear monomode optical fibre, is completely 



integrable [pl| , p% , |33| . The reduced Maxwell-Bloch equations or their gener- 
alisations considered hereafter, as RMB-equations are completely integrable 
too. The 1ST method for both the NLS and the RMB-equations enabling to 
solve some non-linear evolution equations is based on the same spectral prob- 
lem. The model of the USP propagation in a Kerr-type non-linear medium 
doped by resonant impurity atoms incorporates both of these systems. But 
there are no reasons for the resulting system of equations to possess a com- 
plete integrability. 

The evolution of the USP propagating in a non-linear monomode op- 
tical fibre in z- direction is described by the equations, which generalise 



the Maxwell-Bloch equations | 5"9|| . We could name them the non-linear 
Schrodinger and Bloch equations. 

^ + ^ + ^\q\ 2 q + a(p) = 0, (25) 
?l = i5p + 2ifR 3 , ^ = if (q*p - qp*) , (26) 
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where q is the normalised slowly varying complex envelope of the USP defined 
by the following expression 



E(x, y, z, t) = A Q q(z, t)*(ac, y) exp[i((3 z - u t)], 

ty(x, y) is a mode function that determines the transverse distribution of the 
electric field over the fibre cross-section. Here ( = z/L D , t = (t — z/v g )tpQ, 
are normalised independent variables of co-ordinate and time, accordingly, 
t p0 is a pulse duration at z = and v g is the USP propagation group velocity. 
The interaction of the radiation with the resonant impurities is characterised 
by the dimensionless constant / = dAot p o/h, where d is an effective matrix 
element of the dipole transition between the resonant states. The coefficient a 
is expressed in terms of the dispersion length Lp and the resonant absorption 

length L a [Q as a = LdL' 1 f^ 1 , where L a = chn e ff (27ru nAd 2 tpoj 

Now let us consider ultra-short pulses propagation in a Kerr-type dis- 
persive medium when the transition between energy levels of the impurity 
atoms are degenerated over the orientations of the total angular momentum 
ji and 2i- The same system of equations appears when fibre contains three- 
level impurity atoms, so that formally we could speak about optical vector 
solitons in a general case. Let the electric field strength write as 

E (j \x,y,z,t) = Aoqj(z,t)^(x,y) exp[i(P z - u t)). 

The equations for the normalised envelope and the variables of the atomic 
resonant system can be written in a unified form as it was done above 



:% , , ..1*2, 

d( <9r 2 



dpj a) 

dr 



MPj* = -if (£ qMf - Qj N^ , (27) 



£m« = -if (rfpW - q lP f>) , = if E (g*Pf - q 3 P^) 

j 

where j and I mark the spherical components of the vectors, and Pj = 

J2a,/3 a Pj ■ The variables into the Bloch equations in the system ( p7]) have 
been determined above. 
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It has been found [0S| that both system of equations 



and (0) 



are 

integrable ones only under condition L^Lj^ = 2f 2 . This condition implies 
that the soliton of the self-induced transparency should simultaneously be 
also a soliton of the NLS equation. To put this another way, the amplitude 
and duration of the 27r-pulse should precisely be of such values that the 
corresponding self-action (due to the high-frequency Kerr effect) would lead 
to complete compensation of the dispersion broadening of the USP. Thereby, 
the existence of the optical solitons in a fibre doped with resonant impurities 
is restricted. 

Yet another example of integrable model describing USP propagation in 
fibre was proposed in IpOfl . There the NLS equation fl25|) has been replaced 



by the Hirota-like equation: 



.dq ld 2 q 2 (d 3 q 2 dq 3 d\q\ 2 \ 

The resulting system of equation admits the Painleve property for certain 
relation between the physical parameters involved in the model. Also the 
zero-curvature representation of this system was explicitly found. 



4 Conclusion 

In conclusion it will be considered some problems related with an integrable 
models and developments of the ideas mentioned above. 

In the traditional scheme of the 1ST method the spectral parameter of 
the auxiliary linear problem is considered as a constant. The generation of 
the 1ST method, where the spectral parameter depends on time and space 
co-ordinate was proposed in |nj. Non-linear equations arising in this ap- 
proach include the explicit dependence on co-ordinates. In framework of 
this method the generation of the Maxwell-Bloch equations was constructed. 
This new system describes evolution of the USP in two-level medium wherein 
the pumping of the excited states of atoms operated continuously. 

Propagation of the USP in a long two-level amplifier was considered in 
|62| . |63j . The 1ST method was applied to the SIT equations in order to obtain 
its uniform asymptotic solution at long distance. It was shown that the 
amplified pulse is always of a quasi self-similar nature. In the neighbourhood 
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of the wave front USP is described by the Painleve equation, whereas far from 
the front the solution goes into the rapidly oscillating self-similar regime. 

The SIT equations has been employed to describe the phenomenon 
of superfluorescence (i.e., short pulse generation from initial fluctuations of 
polarisation in a high inverted of the population in two-level medium). Con- 
trary to the self-induced transparency, the superfluorescence pulse is gen- 
erating from an unstable state and 1ST method is to be re-formulated for 



boundary problem at the half t-axis. It was done in [64 



Progress in laser physics has made possible the generation of intense ultra- 
short pulses to produce the multi-photon processes. The simplest one devel- 
ops when the resonant matter interacts with a pair of USPs of different carrier 
frequencies uj\ and co> 2 so that uii± u 2 coincides with a transition frequency 
co>2i- The non-linear processes are two-photon absorption and Raman scat- 
tering, respectively. The coherent pulse propagation takes place under both 
resonance conditions when the pulse duration is much less then relaxation 
times. This process is similar to SIT and it was referred to as two-photon 
self-induced transparency (TPSIT). This phenomenon can be described in 
framework of the generalised Maxwell-Bloch equations |65, These equa- 
tions were converted into new system of equation for new variables, which are 
quadratic functions of initial slowly varying envelopes of the electromagnetic 
fields |)B], j57|. As was shown by Kaup and Steudel |B7|, |BS[ under certain 



conditions these equations can be solved by using modified 1ST method. 

It is interesting that both McCall-Hahn equations in the sharp-line limit 
and exact resonance and GSIT equations can be mapped into the model of 
principal chiral fields. The Kaup-Steudel equations are related with the same 
model Ipm . Recently the authors of |70| have investigated the GSIT equa- 
tions and found the hidden non-Abelian group structure of these equations 
in the case of multi-level resonant medium. They have discovered that a 
non-degenerate two-level system of self-induced transparency is associated 
with symmetric space G/H = SU(2)/U(1) while three-level V- or A-systems 
are associated with G/H = SU(3)/U(2). The same symmetric space is 
associated with the degenerate two-level system of SIT in the case of tran- 
sitions jx = 1 — > J2 — and ji = — > ji — 1. When one considers the 
transition j\ — 1/2 — > — 1/2; the GSIT equations are associated with 
G/H = (SU(3)/U(2)) 2 . There are many complex aspects related to the de- 
generation of the energy levels in a three-level system of SIT. For instance, 
the transitions between states j a = j c = 0,jb = 1 (or j a — j c — 1, j b — 
) are associated with the symmetric space G/H = SU(4) / S(U(2) x 17(2)) 
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(accordingly G/H = SU(5)/U(4) ). 

The degenerate four-wave interaction consisting of two counter-propagating 
pulses in a cubic non-linear medium, with arbitrary polarisations which can 
vary through each pulse has been considered in |7J], [72| . The group velocity 
dispersion was neglected that results in the intensity envelopes propagating 
unchanged from their initial forms. However, the polarisation evolves ac- 
cording to the non-linear interaction and it is found it can exhibit soliton 
behaviour. The problem is best approached if one introduce the following 
vectors (i.e., Stokes vectors) 

S (±) = {si^gf ) + c.c, igpgM + c.c, |gW| 2 - \tf f} 

where q^y are the transverse normalised slowly varying envelopes of the 
pulses propagating in the +z and — z directions, respectively. By using the 
cone co-ordinate 2x± = ^p(z ± ct), the evolution equations for can be re- 
written as the equation of motion for chiral field on the group 0(3) : 

«9 + S (+) = S (+) x JS ( - } , <9_S ( - } = -S<"> x JS (+) (29) 

where J diagonal matrix is defined by constants of the model. These equa- 
tions are integrable by 1ST method [173, WM. It should be pointed that in 



the anisotropic case (i.e., matrix elements of are different) the equations (p9|) 
describe the domain wall (DW) in the configuration of the Stokes vectors. As 
far as we know it is the first example of domain structure in non-linear optics, 
where DW's represent space regions with different stable polarisation states 
of the interacting optical waves. The Hamiltonian achieves the minimum 
on these optical domains, and DW's are the regions of polarisation switch- 
ing. Recently, study of the DW's was considered in systems of coupled NLS 



equations governing propagation of light in non-linear fibres ||75|| . An impor- 



tant point is that group- velocity dispersion was taken into account. The two 
types of the DW were found: the one between different elliptic polarisations 
in the bimodal fibre and a dark soliton in one core of the dual-core coupler. 
However, both the model considered in [7JJ and models quoted there are not 
integrable ones. 
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